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High-dimensional biphoton states are promising resources for quantum 
applications, ranging from high-dimensional quantum communications 


to quantum imaging. A pivotal task is fully characterizing these states, 
which is generally time-consuming and not scalable when projective 
measurement approaches are adopted; however, new advances in 
coincidence imaging technologies allow for overcoming these limitations 
by parallelizing multiple measurements. Here we introduce biphoton 
digital holography, in analogy to off-axis digital holography, where 
coincidence imaging of the superposition of an unknown state witha 
reference state is used to perform quantum state tomography. We apply 
this approach to single photons emitted by spontaneous parametric 
down-conversion in a nonlinear crystal when the pump photons possess 
various quantum states. The proposed reconstruction technique allows 
for amore efficient (three orders of magnitude faster) and reliable (an 
average fidelity of 87%) characterization of states in arbitrary spatial 
modes bases, compared with previously performed experiments. 
Multiphoton digital holography may pave the route toward efficient and 
accurate computational ghost imaging and high-dimensional quantum 
information processing. 


Photonic qudits are emerging as an essential resource for environment- 
resilient quantum key distribution, quantum simulation and quantum 
imaging and metrology’. The availability of unbounded photonic 
degrees of freedom, such as time-bins, temporal modes, orbital 
angular momentum (OAM) and radial number’, allows for encod- 
ing large amounts of information in fewer photons than would be 
required by qubit-based protocols (for example, when using only 
polarization). At the same time, the large dimensionality of these 
states, such as those emerging from the generation of photon pairs, 
poses an intriguing challenge for what concerns their measurement. 
The number of projective measurements necessary for a full-state 
tomography scales quadratically with the dimensionality of the Hil- 
bert space under consideration’. This issue can be tackled with adap- 
tive tomographic approaches’ ° or compressive techniques”, which 


are, however, constrained by a priori hypotheses on the quantum state 
under study. Moreover, quantum state tomography via projective 
measurement becomes challenging when the dimension of the quan- 
tum state is not a power of a prime number*. Here we try to tackle the 
tomographic challenge, in the specific contest of spatially correlated 
biphoton states, looking for aninterferometric approach inspired by 
digital holography’, familiar in classical optics. We show that the 
coincidence imaging of the superposition of two biphoton states, 
one unknown and one used as a reference state, allows retrieving the 
spatial distribution of phase and amplitude of the unknown biphoton 
wavefunction. Coincidence imaging can be achieved with modern 
electron-multiplying charged coupled device cameras”, single pho- 
ton avalanche diode arrays” “ or time-stamping cameras’”"*. These 
technologies are commonly exploited in quantum imaging, such as 
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Fig. 1| Biphoton state holographic reconstruction. Pictorial representations 
of equation (2) (inthe simplified scenario of a two-dimensional space (X, X,)) for 
different scenarios in which the biphoton wavefunctions have variable spatial 
correlations. To mimic typical SPDC states, we modelled the reference biphoton 
as a product exp(—a(X; — XI exp(—(Xj + XF I + ik(X; + X,)) andthe 
unknown wavefunction as exp(—a(X; — XYI exp(—(Xj + XF 10) 

h((Xi + Xs)/02), where h,(x) are Hermite polynomials. The parameter a 
quantifies the narrowness of the diagonal correlations. Alongside the three- 


dimensional plot of C(X;, X,), we show the marginal correlation f C(X;,X,)dX; 
(which corresponds to the coincidence image obtained when no spatial 
post-selection is performed) and the section C(X;, X;), which can be obtained 
post-selecting on diagonal correlations. a-c, We chose o = 1, 0, = 0.6, k= 211/ 
(0.20), and a =1 (a), a=100 (b) and a= 400 (c). The correlation width ya is 
reported in each pane. We see that, in the strong correlation limit (c), interference 
is also retrieved in the marginal distribution. 


ghost imaging experiments” or quantum super-resolution”°”, as well 
as for fundamental applications, including characterizing two-photon 
correlations?” , imaging of high-dimensional Hong-Ou-Mandel 
interference”, and visualization of the violation of Bell inequali- 
ties”. Holography techniques have been recently proposed in the 
context of quantum imaging” ”; demonstrating the phase-shifting 
digital holography in a coincidence imaging regime using polari- 
zation entanglement”, and exploiting induced coherence, that is, 
the reconstruction of phase objects through digital holography of 
undetected photons”. 

In this work, we focus on the specific problem of reconstructing 
the quantum state (in the transverse coordinate basis) of two photons 
emerging from degenerate spontaneous parametric down-conversion 
(SPDC). These states are characterized by strong correlations in the 
transverse position (considered on the plane where the two-photon 
generation happens), which can be observed in other kinds of pho- 
ton sources such as cold atoms”. In these sources, the two-photon 
wavefunction strongly depends onthe shape of the pump laser used 
to induce the down-conversion process”. The most commonly used 
approach in the literature to reconstruct the biphoton state emitted 
by a nonlinear crystal is based on projective techniques”. This 
method has drawbacks concerning measurement times (as it needs 
successive measurements on non-orthogonal bases) and the signal 
loss due to diffraction. We proposed an imaging-based procedure 
capable of overcoming both of the issues mentioned above, while 
giving the full-state reconstruction of the unknown state. The core 
idea lies in assuming the SPDC state induced by a plane wave as known, 
and in superimposing this state with the unknown biphoton state. 
Unless the superposition is achieved directly on the crystal plane, 
a full analysis of the four-dimensional distribution of coincidences 
is necessary to retrieve the interference between the two wavefunc- 
tions. This information can be visualized by observing coincidence 
images, defined as marginals of the coincidence distribution obtained 
integrating over the coordinates of one of the two photons. In fact, 
obtaining coincidence images after post-selecting specific spatial 
correlations allows retrieval of the phase information, likewise in 
cases in which the state does not exhibit sharp spatial correlations. 
We demonstrate this technique for pump beams in different spatial 
modes, including Laguerre-Gaussian (LG) and Hermite-Gaussian 
(HG) modes. We investigate several physical effects from the 


reconstructed states, such as OAM conservation, the generation of 
high-dimensional Bell states, parity conservation and radial correla- 
tions. Remarkably, we show how, from a simple measurement, one 
can retrieve information about two-photon states in arbitrary spatial 
mode bases without the efficiency and alignment issues that affect 
previously implemented projective characterization techniques. 
Depending on the source brightness and the required number of 
detection events, the measurement time can be of the order of tens of 
seconds, whereas the previously implemented projective techniques 
required several hours and were limited to the exploration of a small 
subspace of spatial modes. As a latter example, we give a proof of 
principle demonstration of the use of this technique for quantum 
imaging applications. 


Theoretical background 

We start by considering a general scenario in which the superposition 
between two biphoton states is created. We label each state as |,) 
and |w,,), where the subscripts rand u stand for reference and unknown, 
respectively. The state |,) is considered as known (for example, 
it can emerge froma source that was previously characterized), whereas 
the goal is finding |,,) from measurements performed on the super- 
position state |Wyo7) = |r} + |W, To simplify the analysis, we consider 
the case in which the photons are frequency-degenerate and with the 
same polarization. The biphoton states can be thus decomposed in 
transverse spatial degrees of freedom (or, equivalently, in the transverse 
momentum) as 


Wn) = | d?X; d?X, Wy(Xi, Xs) [X;) ® IXs), (1) 


where ņ =r, u; , isa complex function of the transverse coordinates 
ofidler and signal photons (X,, X,) ata given plane. The superposition 


state is 
ror = ff ex ex 


IPX, Xs) + Pu(Xi, Xs)] Xi) @ Xs). 


If y, is known, one can retrieve information about the phase of Y, 
by coincidence measurements in the transverse position basis. 
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Fig. 2 | Experimental set-up and state reconstruction. a, Sketch of the 
experimental set-up: a 405 nm laser in a Gaussian mode (E,ef) enters a Michelson 
interferometer, where an ultraviolet spatial light modulator (UV-SLM) in one arm 
is used to shape and generate the unknown pump field (€p). The interferometer’s 
output is the superposition of the reference and unknown pump field, which is 
then shined ona 0.5-mm-thick type-I BBO crystal. Photon pairs are consequently 
generated and, after being separated through the beam splitter, sent on single 
photon sensor arrays. The experiment was conducted with one camera and the 
figure is just for illustration, a more detailed description and figure of the set-up 
can be found in Supplementary Information. b, By placing the camera in the far 
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field of the crystal and pumping with a large Gaussian beam, we can reconstruct 
the phase-matching function A sinc?’ (ax - ) by direct imaging. The figure shows 
a scan of the phase-matching function with a nonlinear fit yielding A = 93 +2 
counts, a= (9.1 + 0.2) x 10% mm? and Z= 0.30 + 0.02. Error bars are s.d. extracted 
from the Poissonian counting statistics and assuming as mean values the the 
registered counts (reported as data points). c, Experimental correlations in the 
x- and y-coordinates obtained by placing the sensors in the image plane of the 
crystal. d, Example of reconstructed phase and amplitude of a biphoton state 
(represented in inverted HSV colours) when pumping the crystal witha 
superposition of LG modes: LG, + LG,-. 


The coincidence count rate corresponding to a simultaneous detec- 
tion ofan idler photon in X; and a signal photon in X, is proportional to 


CX X2) := |X, X2|Pror) I” 
IW-(X1,X2) + Pu Xa, X2)? 


1WXir XÐ? + Pu X, Xo)? 
+ [W Xi, XDPu(Xı, Xo) + c.c]. 


(2) 


The last equality displays an interference term containing the 
phase difference between the reference and the unknown biphoton 
wavefunctions. 

We note that equation (2) represents an interference pattern in 
the four-dimensional space (X,, X,). Figure 1a gives a pictorial repre- 
sentation of equation (2) (here considering only one coordinate per 
photon) for a general case in which the reference and unknown wave- 
functions have a large spatial correlation spread. Although registering 
€(X;,X;) is enough to obtain phase information, a further visualiza- 
tion of the experimental interference pattern can be obtained by look- 
ing at filtered coincidence images. A coincidence image corresponds 
tothe marginal distribution / C(X;, X,)dX,(a=iors) where the interfer- 
ence term is generically washed out. Interference in coincidence images 
can generally be retrieved by looking at sections of the correlation 
pattern, for example, extracting coincidence images from the quantity 
C(X;,X,)6(X, + X;), where the 6 function represents a narrow filter 
applied on the four-dimensional correlation pattern (Fig. 1 shows the 
case X, =X, given by the red plot). As we will demonstrate in the follow- 
ing section, this operation can be performed by analysing the measure- 
ments of atime-stamping camera. As shown below for the case of SPDC 
filtered coincidence images, besides being useful for data visualization, 
can allow to isolate different contributions of the biphoton state. 


For the interference pattern to have a good contrast, the ampli- 
tudes of w, and w, have to be of a similar order. For instance, the 
phase of an unknown state with strong position correlations would 
be well resolved if the reference has the same spread of the spatial 
correlations, whereas the information hidden in the interference 
term would be more difficult to retrieve if using a reference state 
that is spatially uncorrelated or anticorrelated. All of these forms 
of correlations can be observed, in different propagation planes, 
within the state created in type-I SPDC by a pump beam that is 
well approximated by a plane wave shining a thin crystal”. In fact, 
such a state exhibits sharp correlations in the near field (that is, 
the image plane of the crystal) and sharp anticorrelations in the 
far field, while in intermediate propagation planes, one observes 
wider correlation patterns. Thus, one can, in principle, use this as 
a reference state for measuring any arbitrary two-photon state. 
The practical challenge would be then to find the right proto- 
col for creating the superposition between the unknown and the 
reference state. 

In this work we consider the case in which the unknown state is 
also generated in an SPDC process. In this scenario, ifthe pump beams 
used for inducing the two SPDC states are in phase, the superposition 
can be generated either by mixing the biphoton superposition ona 
beam splitter (see Supplementary Section 3) or directly inducing the 
two SPDC processes in the same crystal. 

The simplest case of the presented scheme arises when both 
the reference and unknown states exhibit sharp position correla- 
tions, as is observed in the case of type-I degenerate SPDC from thin 
crystals. Figure 1b,c shows how the transition toa sharp correlation 
regime allows one to observe interference in the marginal distribu- 
tion without needing post-selection on the spatial correlations. In 
this widely studied limit, which neglects propagation effects in the 
crystal volume, the two photons are created in the same transverse 
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Fig. 3 | SPDC state reconstruction at an intermediate plane. a-d, When moving 
the camera10 cm away from the crystal image plane, the spatial correlations are 
broader (as shown ina, for the case of x-correlations). The image obtained after 
post-selecting only on temporal coincidence (shown in the inset in a) does not 
exhibit any interference as in the crystal image plane. The red and cyan bands 

in aindicate the two-pixel-wide regions selected to analyse coincidence images 
obtained after selecting spatially correlated (b) and anticorrelated (d) photons, 


c 
5.5 
1 
= Oo T 
E 
E äi 
=- | S9 
n 
1 E -rn o 
(0) 
ie) 5.5 o 5.5 


Ce 


i =X, (mm) x= 


i =X, (mm) 


respectively. The correlation image (b) shows again the interference between 
unknown and reference pump fields, allowing us to reconstruct the unknown 
field with off-axis holography (c). The anti-correlation image in d displays the 
characteristic cone shape of the phase-matching function. This function is not 
identical to the one obtained in Fig. 2b, but can be obtained by propagating the 
amplitude of the latter from the far field to the intermediate plane. 


position; this means that the state can be written as (see ref. 31 and 
Supplementary Section 2) 


1») =. f dpe @ |p). 6) 


in the transverse coordinates of the crystal’s image plane, where €,(p) 
isthe pump field onthe crystal plane (the subscript p stands for pump) 
and NW is a normalization factor. The spatial correlation properties 
of biphoton states in this approximation are of strong interest for 
applications in high-dimensional quantum entanglement and quan- 
tum imaging”***°, 

Ifthe thin-crystal approximation holds for both the reference and 
unknown states, equation (2) contains relevant contributions only for 
the case @(X,, X;), with X; = p. These diagonal contributions of the 
coincidence count rate are given by 

C(P.P) = |Ep(P) + Erer) (4) 

Here, Eef is the pump shape used to generate the reference 
SPDC state. Ideally, €,.¢ can be a plane wave or, in practice, a Gaussian 
beam with a large waist. By controlling the reference pump beam 
one can map any interferometric technique that is used in classical 
optics for amplitude and phase reconstruction to the two-photon 
case. In this work, we experimentally implemented off-axis digital 
holography, where the reference beam is a Gaussian beam with a 
tilted wavefront. In off-axis digital holography, with E,er(x,y) = 
Aexp(—(x? + y?)/w2) exp(i2m(x + y)/A), where A is the amplitude of 
the reference beam, w, the waist and 211/A the magnitude of the 
average transverse wave vector components, one has 


2 2 
[Ere + £pl = lEretl? T |Epl 


(x+y) 


wp (ee? + če), 


where c.c. indicates the complex conjugate. From a spatial Fourier 
transform one can hence isolate the term proportional to €, and recon- 
struct the amplitude and phase of the unknown field. The proposed 
scheme can be implemented in two measurement steps: first, the cor- 
relations inthe crystal image plane are measured to confirm the valid- 
ity of the thin-crystal approximation, second, the coincidences 
corresponding to equation (4) are evaluated and the biphoton state 
extracted from the resulting interference pattern. 

Beyond the thin-crystal approximation, one has to reconstruct 
also the contribution of the phase-matching function, whichis in gen- 
eral a function of the form @(X; - X,). The phase-matching contribution 


can be retrieved either by analysing the far field SPDC intensity distri- 
bution or, more rigorously, by post-selecting anti-diagonal correla- 
tions, that is, analysing the coincidence image @(X, —X). 


Experimental set-up and results 

Following the theoretical description of the previous section, we 
experimentally implemented a platform in which, through off-axis 
digital holography, the biphoton state, emitted via SPDC by a type-I 
B-barium borate (BBO) crystal, is reconstructed. In this proof of prin- 
ciple experiment, we generate the unknown and reference SPDC states 
inthe same crystal. A visual scheme of the set-up is reported in Fig. 2a 
(see Methods and Supplementary Section 5 for details). We built a 
Michelson interferometer which allowed the creation of apump beam 
inthe mode €, + Eyer, where the reference mode is a wide Gaussian with 
a tilted wavefront Eref = exp(—?/w?) exp(i2m(x + y)/A), and Epis gen- 
erated with a spatial light modulator (SLM) placed in one of the arms 
of the interferometer. The mean transverse momentum 2v2m11/A is cho- 
sen to maximize the spatial resolution of the reconstructed field, and 
w,is chosen to be larger than the characteristic waist parameter of £p, 
denoted as w,. The interferometer’s output is sent through the BBO 
crystal, and the state of the two photons is recorded using a 
time-stamping camera (Tpx3Cam). The camera comprises a matrix of 
256 x 256 time-stamping pixels of 55 um size and with ~1 ns 
time-resolution. We collected data in the crystal’s image plane, which 
were used to reconstruct the biphoton state in the thin-crystal approx- 
imation. The emitted signal and idler photons were separated by a 
beam splitter and sent into different regions of the camera sensor, 
allowing one to check for coincidences between different pixels. In 
Fig. 2a we represent, for simplicity, the two sensor areas as two inde- 
pendent time-stamping cameras. We could verify the correctness of 
the thin-crystal approximation by observing the spatial correlations. 
This is shown by the sharp, ~1-pixel-wide, spatial correlations observed 
in all of the cases under analysis (see Fig. 2c for an example). The data 
was collected in 1 min of exposure for each spatial mode under analy- 
sis. In particular, we collected both the interference pattern between 
the two states and a coincidence image of down-converted light 
induced by £, only. The former was used to retrieve the phase of the 
state under analysis, whereas the latter already gives the amplitude of 
the biphoton field. By exploiting this reconstruction, we were able to 
fully characterize the biphoton state. An example of the reconstructed 
phase and amplitude of the biphoton state is in Fig. 2d. Moreover, we 
also characterize the amplitude of the phase-matching function. To 
do so, we collect an image of the far field by placing an additional lens 
in front of the Tpx3Cam camera and entering the crystal with a wide 
Gaussian beam. Fitting the collected data with the function sinc*(ax* - J), 
aand (being fit parameters, we obtained a = (9.1 + 0.2) x 10°° mm? and 
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Fig. 4| OAM and radial correlations in SPDC biphotons. Reconstructed field of 
the biphoton state for different OAM-carrying pump beams in LG, =0,lp modes. 
a, The amplitude and phase of the state for different values of the pump OAM £,. 
b, The OAM correlations density plots of generated SPDC photons. It can be seen 
how increasing the pump OAM, the sum of the OAM values for the idler and 
signal photons shift in agreement with the conservation law of equation (6). 


-5 
a” 0 
-5 0 5 


Ps 


Pi 


Pi 


i fi 
d a : igus P = 
2,=0,p,=0 2,=1, Pp =0 2,=2,p,=0 
F= (98.5 + 0.1)% F= (98.90 + 0.06)% F= (98.1 + 0.1)% 
0 0 
g g 
10 10 
o 10 
Pi Pi Pi 
£,=0,p,=1 €,=1, Pp =1 £,=2,P,=1 
F = (93.6 + 0.5)% F = (91.8 + 0.3)% F = (90.0 + 0.2)% 
o Oo o 
2 g g 
10 10 10 
0 p; 10 
€,=1,P,=2 £,= 2, Py= 2 
F = (89.5 + 0.4)% F = (81.9 + 0.4)% F = (81.6 + 0.5)% 
o [0] ie) 
e g g 
10 10 10 
o 10 o 10 Oo 10 
Pi Pi Pi 
£, =0,p,=3 2, =1 P =3 t, =2 p, =3 
F = (82.6 + 0.9)% F = (76.7 + 0.3)% F = (711.2 + 0.3)% 
Oo te) Oo 
g g 


ie) 10 
Pi 


c, Reconstructed biphoton fields obtained by pumping the crystal with LG modes 
changing the ? and p indexes. d, Shows the correlations in the radial number p. In 
all of the plots, the OAM of signal and idler has been fixed to £; = 0, £, = £,. The 
fidelities are obtained assuming the theoretical state calculated in the thin- 
crystal approximation. The error analysis is reported in Methods. 


(= 0.30 + 0.02. From a, we retrieve a value of the crystal length 
Lexp = 40,a/C = 0.56 + 0.01 mm (where c is the speed of light in the 
medium), which is in very good agreement with the nominal value 
Lom = 0.5 mm. The phase-matching fit is shown in Fig. 2b. This extra 


measurement allows one, in principle, to evaluate corrections to the 
quantum state beyond the thin-crystal approximation. We took a meas- 
urement moving the camera away from the image plane of the crystal 
to show that our approach is not necessarily limited to sharply 
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Fig. 5 | Hermite-Gauss correlations. When pumping the crystal with 

Ep, y) = HGm,n (x, yY), we observe biphoton correlations in the basis of HG 
modes that highlight the parity conservation of the SPDC process. Insets show 
the reconstructed biphoton fields from which the correlations have been 


10 ie} 10 


n =m; 


extracted. Upon each plot, the fidelity between the retrieved field and the 
theoretical one in the thin-crystal approximation is reported. The error analysis 
is discussed in Methods. 


correlated photon states. In this plane, the transverse spatial correla- 
tions are broader, and the two-photon wavefunction takes the form 
(|X; — X,| /2)E((X; + X,)/2)(see Supplementary Section 2 for details), 
where @ is related to the phase-matching and € corresponds to the 
pump field propagated from the crystal to the measurement plane. As 
discussed in the previous section and illustrated in Fig. 1a, in this sce- 
nario, the imaging of the singles or coincidences, without any 
post-selection, does not show any interference (this is illustrated in 
inset of Fig. 3a); however, when we look at the coincidence image after 
selecting diagonal spatial correlations, we obtain the marginal 
S dX, CX, X,)6(X; — Xj) = |E(X)I716(0)/? . In our case, € = Ep +E rer 
and we therefore see the interference between reference and unknown 
pump beam (see Fig. 3b). From this pattern, one can reconstruct the 
pump field contribution to the SPDC state, as shown in Fig. 3c. 
This reconstruction yields only the part of the biphoton state associ- 
ated with the pump field; however, the phase-matching can also be 
reconstructed from the same measurement by extracting 
SEX, X,JO(X; + X,)dX, = |£(0)|"|@(X)|?_ (see Fig. 3b). We note that 
in our set-up, the phase-matching functions of reference and 
unknown state are identical, so we only have access to the absolute 
value of this contribution. It must be noted, however, that a set-up in 
which the two SPDC states are created from two independent sources 
would also allow access to eventual phase structures in the 
phase-matching term (see Supplementary Section 3). Alternatively, 
one can extract the phase-matching contribution at different planes 
and infer its phase using the Gerchberg-Saxton algorithm”. We will 
investigate this non-interferometric approach in a future work. Once 
the biphoton state is given, one can extract any desired information 
about this state, for example, correlations in different degrees of free- 
dom, entanglement and the decomposition in arbitrary sets of 
spatial modes. 


Reconstruction of spatial mode correlations 

One of the degrees of freedom of light that has been extensively 
studied for high-dimensional quantum applications, is the OAM. 
Modes in OAM eigenstates are represented by wavefunctions pos- 
sessing a phase term of the form exp(ié@)in position representation. 
Here, gis the azimuthal angle in cylindrical coordinates, and £ € Z 
is the OAM value (along the propagation direction) in units of Å car- 
ried by a photon in sucha state. Given a set of OAM-carrying modes 
(r,@\p, £) :=f,(r) exp(ié@), where f,(r) denotes an orthonormal set 
of radial functions, the probability of detecting the idler photon with 
OAM £, and the signal with OAM £, is (for fixed p; and p,) 


2 
Pie. = | I ERr, O) fp, OSa r), e+?) rdrdg) . (5) 


In particular, when the pump beam carries an OAM equal to , one has 
the OAM conservation law: 


lp = fi + bs, (6) 


which can be immediately deduced from equation (5) and was first 
demonstrated in ref. 32. We investigated this relationship for several 
OAM values by entering the crystal with LG modes*’. These area set of 
modes defined as: 


LG, elr, P) :=(r, Ip, £) œ (r/w)!'L (272 /w?) exp(—(r/w)”) explieg), 


where Lx) are associated Laguerre polynomials. LG modes are cylin- 
drically symmetric modes carrying OAM and with minimal divergence 
in free space”. We analysed the case in which we entered the crystal 
with states having azimuthal index £ € {0, 1, 2, 3, 4}, the results are 
reported in Fig. 4. We show the OAM detection probabilities without 
specifying the radial function, calculated as Pe, e, = = Dp, wee a where 
we restricted the sums over the radial indexes for Pis =0, ... , 10. It is 
evident how increasing the OAM carried by the pump causes the OAM 
correlations to shift in agreement with equation (6) (as also observed 
in, for example, ref. 32). It has been observed that the SPDC state inthe 
OAM basis violates high-dimensional Bell inequalities**”°. In Fig. 4b, 
wealso report the trace distances between the theoretically calculated 


ye and reconstructed (P; 2.) probability distributions, defined as 


D = Yo, o, |Pewe. — P o |/2. The values of D increase with a higher £, 


mainly due to imperfections inthe pump preparations (hence the pump 
beamis better described as a superposition of OAM states). The distri- 
butions P% o are shown in Supplementary Fig. 4. 

The radial index p € Nof LG modes corresponds to the number 
of radial zeros and can be treated as a quantum number*™*., In Fig. 4c, 
experimental results of biphoton states for pump beams, prepared as 
LGmodes, are shown. The coefficients of the SPDC state decomposition 
in LG modes |) = ae : chy fs lPi: ĉi) ® |Ps,@s) were extracted from 
the reconstructed states. Figure 4d shows experimental correlations 
in radial indexes (with OAM indexes of idler and signal fixed as 
€,=0, = €,). When choosing the waist parameter of the decomposition 
to be equal to the pump waist, the correlations are maximized for 
Pis =P» Psi = 0. This canbe understood fromthe similarity between the 
integral expression of oy ‘and the orthogonality relationship of LG 


modes. The fidelity F = | Sa fp cue was evaluated within the 


considered subspaces (where indices p,, are bounded from O to 10 
and £; = 0, £, =£,), where oF represents the measured coefficients 


and cones represents those expected from the thin-crystal 
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Fig. 6 | Image reconstruction. a, Coincidence image of interference between 
areference SPDC state and a state obtained by a pump beam with the shape of 
a Ying and Yang symbol (shown in the inset). The inset scale is the same asin 
the main plot. b, Reconstructed amplitude and phase structure of the image 
imprinted on the unknown pump. 


approximation—both of the coefficients are obtained after normalizing 
the state in the reduced subspace. 

In this approximation, the parity of the spatial modes is also con- 
served”. This effect can be highlighted by considering pump beamsin 
HG modes”? which, on the crystal plane, read 


HGm,,n, x,y) := (x, ylMp, Np) 
x exp( - (+y?) lw )hm, (x/wp) hn, (Y/wp), 


where h„(x) represents Hermite polynomials of order m. HG modes 
form a complete, orthonormal set with even or odd functions along 
thex- or y-directions; the SPDC state can therefore be decomposed as 
1V) = Erin, Kan, IMi Ni) ® IMs, ns) . This basis has been extensively 


MiNi 

studied®***’ and recently considered for biphoton super-resolution 
measurements**. When studying the SPDC correlations on the basis 
of HG modes (chosen with the same waist parameter w = w, of the 
pump), one has the conservation laws: n, = mod(n; + n;,2) and 
m, = mod(m; + m;,, 2), as can be directly inferred from the parity of 
the integrands appearing in the expression of Kins (see Methods and 
ref. 35 for a detailed proof). Figure 5 shows calculated HG mode cor- 
relations for different states (with reconstructed biphoton amplitude 
and phase shown in the insets). The results show an excellent agreement 
with theory; in particular, the parity conservation is evident from the 
chessboard-like correlation patterns. 

Finally, Fig. 6 demonstrates an example of the potential 
applications of biphoton digital holography. The unknown pump 
beam can carry information about an image or be scattered by a 
three-dimensional object. The information about the scatterer is trans- 
ferred to the SPDC state and can be retrieved through our technique 
(Fig. 6b). We show this in the case of off-axis holography, which can 
present limitations for complex structures due to the limited camera 
resolution. These limitations are not related to our proposal and can 
be improved by employing other approaches, for example, on-axis 
phase-shifting digital holography’. 


Conclusion 

Inthis work we introduced anovel approach for reconstructing the spa- 
tial structure of correlated two photons states. Our proposal exploits 
the coherent superposition of two SPDC states and the possibility of 
imaging the amplitude of this superposition with a time-stamping 
camera. The experimental results showed how, froma single measure- 
ment, it is possible to retrieve, in post-processing, alarge amount of 
information about a two-photon spatial state, including correlations in 
different degrees of freedom, entanglement and spatial mode decom- 
position in arbitrary bases. We focused on the simple case of SPDC 


generated from thin crystals and, for different pump fields, analysed 
OAM and parity conservation, high-dimensional entangled states and 
radial mode correlations. The results show the superiority of this tech- 
nique, compared with projective techniques (for example, the ones 
inrefs. 34,42,45), in the context of benchmarking highly correlated 
quantum states. We observe that, if a projective measurement was 
performed on a 11 x 11 subspace, in line with the ones considered in 
this work for radial modes, several days would be required to accumu- 
late the necessary statistics on 121? projections due to the low count 
rates associated with the lossy techniques used for mode projection. 
In comparison, our approach allows us to obtain the necessary data in 
afew minutes, independently from the dimensionality of the subspace 
to be analysed (the latter is only limited by the camera resolution). We 
also note that, with the time-stamping camera, when spatial correla- 
tions are present, a further background subtraction is possible thus 
improving the signal-to-noise ratio of this kind of measurement. We 
thus achieved a three-order-of-magnitude enhancement on the recon- 
struction time with high fidelities for the biphoton states, obtaining 
an average fidelity equal to 87%. The lowest fidelity values are due to 
imperfect pump preparations or an undesired spatially varying phase 
inthe reference beam and not to intrinsic limitations of the technique. 
Although most of our results are based ona specific kind of two-photon 
state, we also showed how it is possible to generalize to states where 
the correlations are not sharp. A key ingredient is to generate reference 
states with spatial correlations that overlap well with the correlations 
inthe unknownstate. Future investigations will be devoted to the gen- 
eralization of this approach to arbitrary two and multiphoton states, 
realizing experiments where the reference and unknown states are 
generated from distant sources. Moreover, we point out that the same 
technique presented here can be applied to measure biphoton states in 
the time-frequency degrees of freedom, where several interferometric 
approaches have been demonstrated*® “**. Besides the quantum state 
reconstruction, future investigations will be devoted to the generaliza- 
tion of our protocol to imaging experiments. 


Online content 

Any methods, additional references, Nature Portfolio reporting sum- 
maries, source data, extended data, supplementary information, 
acknowledgements, peer review information; details of author contri- 
butions and competing interests; and statements of data and code avail- 
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Methods 

Detailed experimental set-up 

A Gaussian beam with a wavelength of 405 nm is produced through 
the second harmonic generation of a 810 nm pulsed Ti:Sa laser 
(Chameleon Vision II), the latter has a pulse duration of 150 fs and 
a repetition rate of 80 MHz. The beam is magnified to ~1cm beam 
waist and sent in the input of a Michelson interferometer. A reflective 
liquid crystal SLM is placed in one arm of the interferometer. At the 
interferometer output, all of the diffraction orders of the SLM, except 
for the first, are filtered out by a slit placed in the SLM’s Fourier plane. 
The slit allows keeping the reference beam (corresponding to the 
beam going through the interferometer arm without the SLM) with 
a different transverse wavevector, thus allowing to perform off-axis 
digital holography while maintaining good interference stability. The 
phase masks applied on the SLM allows the generation of arbitrary 
optical fields by means of the technique introduced in ref. 49. After 
filtering, both beams (each with the power of the order of 100 mW) 
are collimated and sent on the 0.5-mm-thick type-I BBO crystal for 
SPDC generation. The vertically polarized down-converted light is 
collimated by a lens of focal length f= 25 cm, and split into two cop- 
ies, that is signal and idler photons are sent to two separate paths 
(with 50% probability), by a sequence consisting of a half-wave plate 
rotated by 22.5° and a polarizing beam splitter (effectively working as 
an ordinary beam splitter). The two copies are sent on parallel paths 
(by sending them on another polarizing beam splitter and changing 
the polarizations in such a way as to maximize the intensity on one 
output port). The two beams have a lateral shift, that avoids them 
from being overlapped, and are focused, by means of af=50 cmlens 
in the thin-crystal approximation configuration and by a f= 75 cm 
lens for the reconstruction of the state in an intermediate plane, on 
the TPX3CAM sensor. In front of the sensor, a3 nm bandpass filter is 
applied to ensure the frequency degeneracy of the analysed photons. 
The bandpass filter also ensures a sufficient contrast of the interfer- 
ence fringes. In the far field configuration, used to reconstruct the 
phase-matching function, an additional f= 20 cm lens is placed in 
front of the camera ina confocal configuration. 


Data acquisition and analysis 

The data was acquired by collecting SPDC light on the Tpx3Cam for 
1 min for each dataset. We collected data for scenarios in which both the 
reference and pump beams were sent via the BBO crystal, and cases in 
which only the pump beam was sent via the crystal, after blocking the 
reference arm. The acquired data files report the time stamp at which 
the counts were detected (see refs. 17,24,25,50 for more information). 
In our case, as we shine two copies of the SPDC light on different regions 
onthe camera, we can have counts from these regions detected inthe 
same time window. We considered as coincidences the counts from 
the two regions with a time-stamp difference of 5 ns. We analysed the 
spatial correlations from this set of counts, confirming the validity of 
the thin-crystal approximation. A weak constant background in the 
correlation plot is always present (due to dark counts and background 
light) and can be reduced by removing the counts outside the correla- 
tion region. Coincidence images were obtained by plotting the posi- 
tions of the counts selected as coincidences. 

The resulting coincidence images were analysed using stand- 
ard off-axis digital holography’, as described in the main text. The 
decomposition of the reconstructed states in terms of OAM, HG and 
LG modes has been conducted by direct calculations of the expansion 
coefficients in the respective bases. The errors on the fidelities have 
been obtained by repeating the analysis for different state reconstruc- 
tions, where the original coincidence images were modified pixel by 
pixel by random amounts within the uncertainty, given by the square 
root of the coincidences assuming Poissonian statistics. In the main 
text, we report the average fidelity and the standard deviation over 
twenty different realizations. 


Parity conservation 
The coefficients Kpm in the HG expansion of the SPDC state are given 
by the integral 


Kini, = (Mi nil @ (mg, nsl) |W) = 


w ff eee) ta (2) hn, (2) 
see) ie (>) ra, (= )axay. 


We consider the case in which E(x, y) = F(x) Fy(y) where F, and F, are 
even or odd functions of x and y, respectively (which is the case if the 
pump is ina HG mode). We have Chim = lF, x Ir, where 


lr sa f AE) ron, CAG 


Wp 


with /= m,n for €=x, y, respectively. The product h, (=) hy, (=) is 


even/odd if l; + /, is even/odd. Thus, the integral is zero if the parity of 
Fis different than the parity of l + l, hence the conservation law men- 
tioned in the main text. 
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